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. Abstract. Observations of giant radio halos provide unambiguous evidence for the exis- 

Oh' tence of cosmic ray (CR) electrons and magnetic fields in galaxy clusters. The physical 

mechanism generating radio halos is still heavily debated. We critically discuss the pro- 
^ , posed models for the radio halo emission and highlight the weaknesses underlying each 

explanation. We present an idea how the interplay of CR propagation and turbulent advec- 
tion selects a bimodal spatial CR distribution that is characteristic for the dynamical state 
of a cluster. As a result, strongly turbulent, merging clusters should have a more centrally 
concentrated CR energy density profile with respect to relaxed ones with very subsonic 
turbulence. This translates into a bimodality of the expected diffuse radio and gamma ray 
emission of clusters. Thus, the observed bimodality of cluster radio halos appears to be a 
natural consequence of the interplay of CR transport processes, independent of the model 
. of radio halo formation, be it hadronic interactions of CR protons or re-acceleration of low- 

OO ' energy CR electrons. 

in 

, Key words. Galaxies: clusters: intracluster medium - Astroparticle physics - Gamma rays: 

■ galaxies: clusters - Radio continuum: galaxies - Acceleration of particles - Magnetic fields 

1. Introduction sion, forming the so-called cluster radio halos. 

• 5—1 , Several CRe injection sites can also be identi- 

X ■ Relativistic particle populations, cosmic rays flgj the same synchrotron radiation mech- 

H I (CRs), are expected to permeate the intra- anism: shock waves from structure formation, 

. P. • cluster medium (ICM). Cosmic ray electrons active galactic nuclei (AGN), and winds or 

(CRes) are direcdy visible in many galaxy stripping from cluster galaxies. All these 

clusters via their radio synchrotron emis- should also be injection sites for CR protons 

(CRps) and heavier relativistic nuclei. Due to 

Send ojfprint requests to: C. Pfrommer their higher masses with respect to the elec- 
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trons, protons and nuclei are accelerated more 
efficiently. In our own Galaxy, the ratio of the 
spectral energy flux of CRps to CRes between 
1. . . 10 GeV is about one hundred. Similar ra- 
tios are also expected at least for the injection 
from galaxies and structure formation shock 
waves for the same kinematic reasons.^ 

Cluster CRps should have accumulated 
over cosmic timescales since the bulk of them 
is unable to leave through the persistent infall 
of matter onto the cluster and due to the long 
CRps' radiative lifetimes in the ICM of the or- 
der of an Hubble time throughout the entire 
ICM. CRes suffer much more severe energy 
losses via synchrotron and inverse Compton 
emission at GeV energies, and Bremsstrahlung 
and Coulomb losses below 100 MeV. CRes 
with an energy of ~ 10 GeV emit GHz syn- 
chrotron waves in //G-strength magnetic fields. 
Since the associated inverse Compton and syn- 
chrotron cooling time is Tic,syn ~ 10** yr, these 
CRes must have been recently injected or re- 
accelerated. 

2. Observations of radio halos 

Cluster radio halos are our primary evidence 
for the existence of CRs in galaxy clusters. 
They are spatially extended regions of diffuse 
radio emission, which have regular morpholo- 
gies (resembling the morphology of the X- 
ray emitting thermal ICM plasma). Their radio 
synchrotron emission is unpolarised, due to the 
contribution of various magnetic field orienta- 
tions along the line of sight, and Faraday rota- 
tion de-polarisation. 

Cluster radio halos come in two sizes: clus- 
ter wide and therefore giant radio halos and ra- 
dio mini-halos. The former are predominantly 
found in clusters showing merger activities 
whereas the latter are found in very relaxed 
clusters which developed a cool core that har- 
bors the mini-halo. The radio (mini-)halo lumi- 
nosity correlates with the X-ray emissivity of 
the cluster (see Fig. 1). A large fraction of clus- 
ters do not exhibit significant radio halo emis- 
sion, and only upper limits to their synchrotron 

' For an extended list of references and a more 
detailed discussion, see EnBlin et al. (2011) which 
our proceeding closely follows. 



flux are known. About half of the radio defi- 
cient clusters, for which we have Chandra data, 
show clear evidence for some level of cool core 
structure (Kq < 40 ke V cm^) as can be seen in 
Fig. 1. This could either imply that these clus- 
ters are in the intermediate state between hav- 
ing giant radio halos because of merging activ- 
ity and having mini halos due to strongly devel- 
oped cool cores. On the other hand there could 
be two populations of clusters - cool cores and 
non-cool cores - and the corresponding radio 
luminosity responds sensitively to the level of 
injected turbulence by either AGN or cluster 
mergers, respectively. 

2. 1. Hadronic models 

In the hadronic model the accumulated CRps 
continuously inject radio emitting CRes into 
the ICM due to well known hadronic process 
PCR+P -^7:"^ + ...^ e'^+vJvc + Vf,+Vfj + ... 
The hadronic model has advantages: 

- AH required ingredients are available: am- 
ple sources of CRps (structure formation 
shocks, AGN, galactic winds), gas protons 
as targets, magnetic fields. 

- Smooth and regular morphology of halos 
are a consequence of the long lifetime of 
CRps which implies a volume fiUing clus- 
ter distribution. 

- Using analytical arguments and hydrody- 
namical simulations, the predicted lumi- 
nosities, scaUngs (Ly - Lx), and morpholo- 
gies match observations without tuning 
(Miniati et al. 2001; Pfrommer et al. 2008; 
Kushnir et al. 2009). 

- The model predicts power-law spectra as 
observed. 

There are also issues with the hadronic model: 

- About two thirds of the most X-ray lumi- 
nous clusters do not exhibit radio halos, 
whereas the hadronic model seems to sug- 
gests that all clusters exhibit halos [will be 
addressed in the following]. 

- The model does not explain all reported 
spectral features (curvature, spectral steep- 
ening) [wiU be addressed in the following]. 
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Fig. 1. Correlation of radio halo luminosities with cluster properties (EnBlin et al. 2011). Left: 
Radio halo luminosity vs X-ray luminosity. Right: Radio halo luminosity vs central entropy 
indicator Kq for the subsample of clusters for which high resolution Chandra data are available. 



The hadronic model makes a testable predic- 
tion: the radio halo emission should always 
be accompanied by weak diffuse gamma-ray 
emission, due to the hadronic production of 
neutral pions and their decay into gamma-rays, 
Pgr + P 7r" + ... — » 2y + . . . The current 
upper Umits on diffuse gamma-ray flux from 
cluster of galaxies by the Fermi collaboration 
(Ackermann et al. 2010) are still well above 
the predictions of expected fluxes, even for 
the most optimistic assumptions about the CR 
acceleration efficiency (Pinzke & Pfrommer 
2010). They are far off the minimal gamma- 
ray flux expected in the limit of strong mag- 
netic field strength (» 3/iG; Pfronomer 2008; 
Aleksicetal. 2010). 

2.2. Re-acceleration models 

In re-acceleration models, a pre-existing CRe 
population at lower energies of about 0.1-10 
GeV gets re-accelerated into the radio emit- 
ting regime of about 10 GeV by plasma waves. 
These are generated by the turbulence during 
and after a cluster merger event. Some level 
of re-acceleration has to happen most of the 
time or frequently enough in order to prevent 
the CRe population in the cluster center from 
loosing its energy completely due to Coulomb 
losses on a timescale of about 1 Gyr. The ad- 
vantages of the re-acceleration model are: 



- All required ingredients are available: ra- 
dio galaxies and reUcs to inject CRes and 
plasma waves to re-accelerate them. 

- The bimodality of radio halo luminosities 
is explained by the presence and decay of 
the re-accelerating turbulence in merging 
and relaxed clusters, respectively (Brunetti 
et al. 2009). 

- Reported complex radio spectra in some 
clusters emerge naturally by interplay of of 
acceleration and cooling. 

The issues with the re-acceleration model are: 

- Fermi 11 acceleration is inefficient and 
scales with Wwave/<^^ ^ 1' efficiency in cur- 
rent models is fitted to explain data and not 
derived from first principles. 

- Current models neglect advective energy 
losses by waves that propagate outwards 
and dissipate in the outer regions. 

- Intermittency of turbulence might be diffi- 
cult to reconcile with the observed regular- 
ity of radio halos. 

- Observed power-law spectra require fine 
tuning. 

- CRes cool rapidly in the central regions on 
timescale of 1 Gyr (for He = 3 x 10"^ cm~^) 
[will be addressed in the following]. 

A testable prediction with upcoming sensitive 
radio telescope arrays is that low X-ray lumi- 
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Fig. 2. Sketch of the interplay of CR streaming and turbulent advection for a single flux tube 
in a stratified atmosphere with gravity pointing downwards. Left: The dense CRs at the center 
stream along the tube towards the CR depleted regions at larger atmospheric height. Middle: 
CR streaming stops as soon as a homogeneous CR space density is achieved. A turbulent eddy 
(represented by its angular momentum axis) starts to turn the magnetic structure upside down. 
Right: The former outer parts of the flux tubes are compressed at the center, and harbor now 
an overdense CR population, whereas the former inner parts are expanded at larger atmospheric 
scale height and therefore have now an underdense CR population. Again CR streaming sets in. 



nous clusters should not exhibit radio halos 
(Cassano et al. 2008). 



3. Cosmic ray transport 

3.1. Confined cosmic rays 

To begin, we consider an isolated magnetic 
flux tube with CRs confined to it to illustrate 
the interplay of advection and streaming with 
a basic picture. This represents the limiting 
case of confined CRs and will be generalised 
in the next section. Imagine a magnetic flux 
tube frozen into the plasma which is distributed 
in a stratified pressure atmosphere of a clus- 
ter as shown in Fig. 2 on the left. Any central 
concentration of CRs will escape due to CR 
streaming on a timescale of Tst - Lg/v^t, where 
Lb is the magnetic bending scale and i/jt the 
CR streaming velocity along the magnetic field 
which is of order the sound speed in the clus- 
ter plasma.^ This leads to a homogeneous CR 

^ In a low-yS plasma, the CR streaming velocity 
is linked to the Alfven velocity, which exceeds the 
sound speed there. However, this can obviously not 
be true in a high-yS plasma. It would imply that in 
the limit of vanishing magnetic field strength the 
CRs get completely immobile due to the vanish- 
ing Alfven speed. However, for disappearing mag- 
netic fields, the coupling of CRs to the plasma gets 



distribution within the flux tube (Fig. 2, mid- 
dle). Turbulence turns the magnetic structure 
upside down on half an eddy turnover time. If 
this is comparable to, or less than, the CR es- 
cape time, 

— =j,,~ 0(1), (I) 

a good fraction of the CRs from larger radii 
will be compressed towards the center, from 
where they again start streaming to larger 
radii. The transonic turbulence is therefore 
able to maintain a centrally enhanced CR 
density by pumping expanded CR populations 
downwards. As soon as the turbulent velocities 
become significantly subsonic, this pumping 
becomes inefficient, since the streaming will 
be faster than the advection. At this point a 
nearly constant volume density of CRs estab- 
lishes within a closed flux tube, meaning that 
most CRs are residing at larger cluster radii. 
Depending on the level of turbulence, we ob- 
tain either a CR distribution that is peaked to- 



weaker and therefore the CRs should stream faster. 
Thus, there must be a characteristic velocity, below 
which the Alfven velocity is not limiting the stream- 
ing velocity any more. Plasma physical arguments 
indicate that this is roughly the sound speed (Felice 
& Kulsrud 2001; EnBlin et al. 201 1). 
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Fig. 3. Left: CR density profiles for y,u = 1,3, 10, 30, and 100 (from bottom to top at small 
radii) including the same number of CRs each. Profiles are normalised to £'(0)|y, -co- Also the 
more narrow gas density profile is shown (thick grey fine). Right: CR normafisation profiles for 
the same parameters and the gas density profile for typical cluster conditions. 



wards the center or a homogeneous CR distri- 
bution. 



3.2. Mobile cosmic rays 

In reality, CR diffusion perpendicular to the 
mean magnetic field enables CRs to change 
between magnetic flux tubes and thereby find 
paths to more peripheral regions. The accessi- 
ble distance is determined by the level of tur- 
bulent pumping, magnetic topology, and avail- 
able time to stream. In principle, CRs can even 
reach the outskirts of galaxy clusters, where the 
infall of matter onto the cluster behind the ac- 
cretion shocks prevents further escape which 
motivates our term mobile CRs. 

We assume a power-law CR spectrum. 



f{r,p,t) = C{r,t)p-", 



(2) 



where a ~ 2.1 - 2.5 is the spectral index 
and C(r, t) the spectral normalisation constant. 
First, we derive the equiUbrium profile that 
CRs attain if turbulent advection dominates the 
CR transport (and CR streaming is negligi- 
ble). To this end, we assume that (i) the clus- 
ter is characterised by a mean pressure pro- 
file and (ii) that CR propagation operates on 
small scales, permitting CR exchange between 
nearby gas volume elements, but not on large 
scales. Whenever two volume elements come 



close, CRs can be exchanged which establishes 
a constant CR population in any given radial 
shell. During radial advective transport from 
radius r to r', the ICM gas with the entrained 
CRs is compressed or expanded by a factor 
X{r r') = (P(r')/P(r))i/T, where P(r) is the 
pressure profile and 7 = 5/3. The CR rest-mass 
density g{r) = m fdpf{r,p) thus establishes 
- under the influence of advection alone - a 
profile according to 



gir) 



Pir) 
Po 



(3) 



where j]{r) = (Pir)/Poy'^ is flie advective CR 
target profile. 

The CR continuity equation for g in the ab- 
sence of sources and sinks can be written as 



(4) 



with V = v^d + i^di + ■^st the CR transport veloc- 
ity, which is composed of the advective (Uad), 
diffusive (t^di), and steaming (Vst) transport ve- 
locities. These are defined by 



Vg 

IV £.1' 



-/Cdi Vln(£)), 



(5) 
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where k^h is the macroscopically averaged CR 
diffusion coefficient and the passive, advective 
transport via turbulence can be described by an 
additional diffusion process with diffusion co- 
efficient /ftu = LtuVtu/^- We note that the ap- 
pearance of the target density profile r](r) in the 
gradient for ensures that any deviation of 
the CR distribution from target density causes 
a restoring term towards this equilibrium con- 
figuration. The CR space density becomes sta- 
tionary for t; = 0, and this reads in spherical 
symmetry with radially outstreaming CRs 

Vst = '^tu -J- In I - 1 + Kdi ^ Inijo). (6) 
or \r]l or 

Enfihn et al. (201 1) provide an analytical solu- 
tion of this equation which is shown in Fig. 3 
for different values of ytu- 

4. Radio and gamma-ray bimodality 

The gamma-ray emissivity and luminosity of a 
power law CRp spectrum as in Eqn. (2) is 

Sy oc C Qg^, andLy = J dVey. (7) 

The radio luminosity in the hadronic model is 

Sy ocC^pgas , and Lv = I dVsy. (8) 

SB + Bph J 

As shown in Figs. 4 and 5, Ly and Ly inherit 
the strong dependence on the advective-to- 
streaming-velocity ratio, ytu - Vtu/fst- Thus, a 
rapid drop in radio luminosity after the turbu- 
lent merger phase by one order of magnitude or 
more is actually expected in the hadronic halo 
model on a timescale of 0.1-1 Gyr, depending 
on magnetic topology and the macroscopic CR 
streaming speed (EnBUn et al. 201 1). 

5. Conclusions 

CR streaming (and CR diffusion) aims at es- 
tablishing a spatially flat CR profile; hence ex- 
plaining why radio halos are not found in ev- 
ery cluster. CR advection tends to produce cen- 
trally enhanced CR profiles. Thus, CR advec- 
tion and streaming are counteracting transport 



mechanisms. Whenever the former dominates, 
centrally enhanced profiles are established, and 
whenever streaming is more important, a flat 
profile results. 

During a cluster merger, advective veloc- 
ities in galaxy clusters are comparable to the 
sound speed and drop when the cluster re- 
laxes after the merger. Plasma physical argu- 
ments suggest that the microscopic CR stream- 
ing velocity in clusters might of the order of 
the sound speed. Macroscopically it is reduced 
due to magnetic trapping of CRs in flux tubes 
(which is larger for a stronger turbulence) and 
slow cross field diffusion required to escape. 

As a result of this, merging clusters should 
have a much more centrally concentrated CR 
population ffian relaxed ones. This leads natu- 
rally to a bimodality of their gamma-ray and 
radio synchrotron emissivities due to hadronic 
interactions of CR protons. Also in the re- 
acceleration model of cluster radio halos these 
transport processes should be essential, since 
the re-accelerated CR electron populations in 
the dense cluster centers is probably too vul- 
nerable to Coulomb losses, to survive periods 
without significant re-acceleration. Transport 
of the longer living electrons at the cluster 
outskirts into the cluster center during cluster 
merger would circumvent this problem. 

We also expect an energy dependence of 
the macroscopic CR streaming speed, which 
then should lead to a spatial differentiation of 
the spectral index of the CRp population and 
any secondary radio halo emission. Such spec- 
tral index variation in the radio halo should be- 
come especially strong during phases of out- 
streaming CRps, i.e. when a radio halo dies due 
to the decay of the cluster turbulence. 



References 

Ackermann, M., Ajello, M., AUafort, A., et al. 

2010, ApJ, 717, L71 
Aleksic, J., Antonelli, L. A., Antoranz, R, et al. 

2010, ApJ, 710, 634 
Brunetti, G., Cassano, R., Dolag, K., & Setti, 

G. 2009, A&A, 507, 661 
Cassano, R., Brunetti, G., Venturi, T., et al. 

2008, A&A, 480, 687 



Pfrommer et al.: Cosmic ray transport in galaxy clusters 7 




Fig. 4, Left: Gamma-ray emissivity profiles for the CR distributions in Fig. 3 and X-ray emissivity 
profile of the ICM in grey. Emissivities are normalised to the central emissivity of a cluster with 
ytu = cx). Right: Total gamma-ray flux due to hadronic CRp interactions with the ICM nucleons 
as a function of ytu = ^st/'^'ad and for a = 2.25, 2.5, and 2.75 (soUd, dashed, and dashed-dotted 
lines, respectively). Normalised to Ly for = oo. 
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Fig. 5. Left: Radio emissivity profiles for the cluster shown in Fig. 3 assuming the same mag- 
netic field profiles with Bq - 6/iG and <5b = 0. Emissivities are normalised to the central radio 
emissivity of a cluster with ytu = °°- The X-ray profile is shown in grey. Right: Total radio 
flux due to hadronic CRp interactions with the ICM nucleons as a function of ytu = T"st/Tad and 
for different dependencies of the magnetic energy density on the turbulence level, bb k n(r)y*|f 
and parametrised by Bg = 0, 1, and 2 (solid, dashed, and dashed-dotted Unes, respectively). 
Normahsed to L, for 7^ = 10. (Central field strength Bq = 6juG, a = 2.5). 
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